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Reminder 



• Energy 

■ 1 J = 1 m 2 • kg • s" 2 , the joule (What is llcT?) 

• Potential 

■ 1 V = 1 m 2 • kg • s-3 • A" 1 , the volt (l eV? lfcT/e?) 

■ (note that 1 s • A = 1 C, the coulomb) 

• Force 

■ A vector 

■ 1 N = m • kg • s -2 , the newton 



Potential contours 




What would happen to a positive test charge? 
Is it favourable to move here from infinity? 



Potential contours 




Phenomenological 
model (circa 1785) 
for charge-charge 
interactions in a 
vacuum 

Relates potential to 
charge for 
homogeneous 
dielectric materials 
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('unitless', that is, with 
dimensions of l) 
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Vacuum permittivity 
(8.854 x 10 12 C 2 ] 1 m 1 ) 
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Poisson equation 

• Generalize Coulomb's law to inhomogeneous 
dielectric coefficient 

-£VV(x) = 4/r/(x) >-V-/?(x)V^(x) = 4;r/(x) 

• Describes electrostatic potential due to: 

■ Inhomogeneous dielectric 

■ Charge distribution 

• No mobile ions 

-V-£(x)V0(x) = f(x) 

= J]q i S(x-x i ) 

i 

^(oo) = 0 



Molecular 



Multiple dielectric values: 

■ 1: vacuum 

■ 2 to 4: atomic polarizability 
(solid) 

■ 4 to 10: some liberation, minor 
sidechain rearrangement 

■ 10 to 20: significant internal 
rearrangement 

Multiple surface definitions 

■ Van der Waals 

■ Splines 

■ Molecular surface 



functions 




Poisson equation: energies 



• Total energies obtained from 

■ Integral of polarization energy 

■ Sum of charge-potential interactions 
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Poisson equa 

Total energies obtained 
from 

■ Integral of polarization 
energy 

■ Sum of charge-potential 
interactions 

Numerically-calculated 
energies contain self- 
interaction terms: 

■ Infinite (for analytic 
solution) 

■ Very unstable (for 
numerical solution) 

Self-interactions must be 
removed 



ion: energies 



For Coulomb's law 
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Solvation 

Solvation energies 
obtained directly 
from reaction field 

Difference of 

■ Homogeneous 

■ Inhomogeneous 
dielectric calculations 

Self-energies removed 
in this process 



energy 

\ olv G = G[0 2 ]-G[0 2 ] 

^ i 
^ i 



Screening of electrostatic potential 



Think about mobile ions 

Salt concentration in solution makes the 
charge on the protein 'less effective' 



No mobile ion 




Mobile ion: screening 




Poisson-Boltzmann equation 



• Screening: Debye-Hlickel theory 

• How do the mobile ions organize themselves? 

■ Try the Boltzmann distribution 

• Result: Nonlinear partial differential equation 

-V -a (x) V j (x) - 4^ O^e-"*'^ = q>S(x ~ h ) 

m i 

^(oo) = 0 



Coefficients: charge distribution 



m 



Charges are delta 
functions: hard to model 

Often discretized as splines 
to 'smooth' the problem 




Coefficients: mobile ion distribution 

-V • e (x) V <f> (x) - 4^ Q m n m e-^)^) = 4 ^ ^ ( x _ X; ) 

• Provides: 

■ Bulk ionic strength 

■ Ion accessibility 

• Usually constructed 
based on 'inflated 
van der Waals radii' 




Coefficients: dielectric function 

m i 

Describes change in 
dielectric response: 

■ Low dielectric interior (2 to 20) 

■ High dielectric solvent (80) 

Many definitions: 

■ Molecular (solid line) 

■ Solvent-accessible (dotted line) 

■ van der Waals (gray circles) 

■ Inflated van der Waals 
(previous slide) 

■ Smoothed definitions (spline- 
based and Gaussian) 

Results can be very sensitive 
to the choice of surface!!! 




Poisson-Boltzmann energies 



G M = 7-f[/^-f( v ^)+^ 2 ( cosh M - l ) 

Fixed charge- Dielectric Mobile charge 

potential interactions polarization energy 



dx 



Discretization: local methods 

• Polynomial basis functions (defined on interval) 

• 'Locally supported' on a few grid points 

• Only overlap with nearest-neighbors sparse matrices 




Boundary element Finite element Finite difference 

(Surface discretization) (Volume discretization) (Volume discretization) 



Input and output: blackbox 



• Software: APBS, WHAT IF, etc. 

• Input: 

■ Solute 'pqr': atom positions, charges, and radii 

■ Solvent: dielectric constant, spatial accessibility 

■ Mobile ions: ionic strength, spatial accessibility 

• Output: potential contours, energies, pK a 

• Use thermodynamic cycles to obtain the 
desired energy values 



Solvation energy 

Physical model: transfer solute from dielectric of 1 (vacuum) 
to 80 (water) 

Computational model: transfer solute from homogeneous 
dielectric to inhomogeneous dielectric - eliminate self- 
interaction terms 

Two models can be reconciled through free energy cycle: set 
the 'reference state' 



Bindin 



This calculation assumes no 
conformational change! 

Separate calculation into two 
steps: 

■ Calculate electrostatic 
interaction in 
homogeneous dielectric 
(Coulomb's law) 

■ Calculate solvation energy 
change upon binding 
(Poisson or Poisson- 
Boltzmann equation) 

Self-interactions are removed 
in solvation energy calculation 



energy 




0 = AG 1 + AG 2 -AG 3 - 
AG 3 = ( AG 2 - AG 4 ) + AG l 
= AAG solv + AG couI 



Conformational chan 



Same concepts as binding 
energy calculation: 

■ Calculate electrostatic 
energy due to 
configuration change in 
homogeneous dielectric 
(Coulomb's law) 

■ Calculate electrostatic 
energy due to change in 
solvation between 
configurations (Poisson or 
Poisson-Boltzmann) 



Conformation A 




pK calculations 




• HA + H 2 0 H 3 0 + + A" AG a = -kT In K 



• piC a = -log 10 K 



pK a calculations 

• Want acid dissociation 
constant for residues in 
a particular structural 
context 

• Use 'model' pK a s for 
amino acids 

• Calculate pK a from two 
'binding 1 calculations: 

■ Binding of 
unprotonated residue 

■ Binding of 
protonated residue 



Amino acid 


oc- 

Carboxylic 
acid 


a-Amino 


Side 
chain 


Alanine 


2.35 


9.87 




Areinine 

o 


2.01 


9.04 


12.48 


Asparagine 


2.02 


8.80 




Aspartic Acid 


2.10 


9.82 


3.86 


Cysteine 


2.05 


10.25 


8.00 


Glutamic Acid 


2.10 


9.47 


4.07 


Glutamine 


2.17 


9.13 




Glycine 


2.35 


9.78 




Histidine 


1.77 


9.18 


6.10 


Isoleucine 


2.32 


9.76 




Leucine 


2.33 


9.74 




Lysine 


2.18 


8.95 


10.53 


Methionine 


2.28 


9.21 




Phenylalanine 


2.58 


9.24 




Proline 


2.00 


10.60 




Serine 


2.21 


9.15 




Threonine 


2.09 


9.10 




Tryptophan 


2.38 


9.39 




Tyrosine 


2.20 


9.11 


10.07 


Valine 


2.29 


9.72 
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Hands-on tutorials on the web 



• Practical: 

■ http : / / en . wikiversity .org/ wiki /Pois son- 
Bo It zmann_pr of ile_for_an_ion_channel 

• Adaptive Poisson-Boltzmann Solver (APBS) 

■ http : // apbs . source forge . net/ 

■ Prof. Nathan Andrew Baker 

• pK a calculations (WHAT IF) 

■ http : / /enzyme .ucd.ie/ Science /pKa/ 

■ Dr Jens Erik Nielsen 



